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ABSTRACT 

We investigate systematic uncertainties in determining the profiles of the 
solar sound speed, density, and adiabatic index by helioseismological techniques. 
We find that rms uncertainties — averaged over the sun — of ~ 0.02% — 0.04% 
are contributed to the sound speed profile by each of three sources: 1) the 
choice of assumed reference model, 2) the width of the inversion kernel, and 3) 
the measurement errors. The rms agreement between the standard solar model 
sound speeds and the best helioseismological determinations is about 0.07%. The 
profile of the adiabatic index Fi is determined to an accuracy of about 0.02% 
with the MDI data set. The density profile is about an order of magnitude less 
well determined by the helioseismological measurements. Five state-of-the-art 
models, each with a significant difference in the input physics or a parameter 
choice, all give comparably good agreement with global helioseismological 
measurements. We consider four deficient solar models that are constructed 
using either old input data, assuming the ^He -|- ^He fusion reaction does not 
occur, neglecting element diffusion, or artifically mixing the interior of the 
sun. When used as reference models in the inversion process, these deficient 
models yield sound speeds for the sun that differ only by 0.1% from the sound 
speeds obtained using the standard model. We conclude that even relatively 
crude reference models yield reasonably accurate solar parameters. Although 
acceptable for most purposes as reference models, non-standard solar models in 
which the core is artifically mixed or in which element diffusion is neglected are 
strongly disfavored by the p-mode oscillation data. These non-standard models 
produce sound speed profiles with respect to the sun that are, respectively, 4.5 
and 18 times worse than the agreement obtained with the standard solar model. 
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1. Introduction 

Helioseismology has revolutionized our knowledge of the Sun and enriched, stimulated, 
and largely validated our understanding of the evolution of main sequence stars with masses 
comparable to the sun. The statistical and measurement errors in the best modern samples 
of helioseismological frequencies have been reduced to tiny fractional values, ~ a few times 
10~^, which lead to fractional errors in, e. g., the inferred sound velocities that are formally 
of order 10""^ or less. 

We therefore concentrate in this paper on systematic uncertainties. The main result of 
this paper is that the systematic uncertainties in determining the sound speeds are about 
an order of magnitude larger than the statistical errors. 

We need to know quantitatively the accuracy of the solar sound speed and density 
profiles that are inferred from helioseismology since these characteristics are often used, e. 
g., to test opacity calculations (Korzennik and Ulrich 1989, Tripathy, Basu & Christensen- 
Dalsgaard 1997), to investigate equations of state (Ulrich, 1982; EUiott & Kosovichev 1998; 
Basu, Dappen & Nayfonov 1999), and to derive abundance profiles in the Sun (Antia & 
Chitre 1998; Takata & Shibahashi 1998). Moreover, the precise agreement between the 
inferred sound speeds and those calculated with a standard solar model used to predict solar 
neutrino fluxes is a strong argument that solar neutrino experiments require new physics, 
not revised astrophysics (Bahcall, et al. 1997; Bahcall, Basu, & Pinsonneault 1998). 

Solar oscillation data have been inverted to determine the sound speed and density 
profiles over nearly the entire Sun (cf. Dziembowski, Pamyatnykh & Sienkiewicz 1990; 
Dappen et al. 1991; Antia & Basu 1994; Kosovichev et al. 1997), as well as the 
adiabatic index Ti (cf. Antia & Basu 1994; EUiott 1996; Elliott & Kosovichev 1998). 
However, helioseismological determinations of solar parameters generally proceed through a 
linearization of the equations of stellar oscillations around a theoretical reference model of 
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the sun. We are therefore naturally led to ask the question which is reflected in the title of 
this paper: How dependent are the inferred solar characteristics upon the assumed reference 
model? 

Here is how we have answered this question. We have constructed a broad range of 
conceivable modern solar models, each model with different input physics or assumptions, 
and used this broad set of models as reference models to calculate solar sound speeds and 
densities with different data sets. The range of inferred sound speeds and densities deflne 
empirically the systematic uncertainties that are inherent in using references models to 
invert helioseismological data. 

Five of the solar models are roughly comparable to each other and use physics that is 
within the currently acceptable range; these models are deflned in § |^, where they are given 
the names STD, PMS,ROT,R78,R508. We also use four different models, each of which is 
deflcient in one or more major aspects of the input physics. The deflcient models are also 
deflned in § ^, where they are named OLD, S34 = 0, NODIF, and MIX. 

To explore the systematic uncertainties associated with choosing a particular data set, 
we have chosen three different sources of observational data. We have used solar frequencies 
obtained by: 1) the Michelson Doppler Imager (MDI) instrument on board the Solar and 
Heliospheric Observeratory (SOHO) during the flrst 144 days of its operation (cf. Rhodes 
et al. 1997); 2) the set of frequencies obtained by observations of the Global Oscillation 
Network Group (GONG) between months 4-14 of its observations; 3) a combination of the 
data from the Birmingham Solar Oscillation Network (BiSON; cf. Chaplin et al. 1996) and 
the Low-£ (LOWL; cf. Tomczyk et al. 1995) instrument. The third set is same as the one 
used by Basu et al. (1997) and is described in detail there. The MDI and GONG sets have 
good coverage of intermediate degree modes. The MDI set has p-modes from £ = up to a 
degree of £ = 194 while the GONG set has modes from £ = up to £ = 150. However, both 
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these sets are somewhat deficient in low degree modes. The BiSON+LOWL combination, 
on the other hand, has a better coverage of low degree modes, but has modes from / = 
only up to £ = 99. 

We concentrate in this paper on the properties of helioseismological inversions. We 
discuss the implications for neutrino physics in Bahcall, Pinsonneault & Basu (1999, in 
preparation) and explore the results of mixing and rotation on element abundances and 
helioseismology in Pinsonneault, Basu & Bahcall (2000, in preparation). 

The present paper is organized as follows. We summarize briefly the inversion technique 
in § 1^ and then describe the solar models used in § |^. Before investigating the systematic 
uncertainties due to the choice of reference model, we first investigate in § ^ the uncertainties 
due to the choice of data set and in § ^ the uncertainties due to the finite resolution of 
the inversion kernel. We present the inferred solar profiles of sound speed and density in 
§ 1^ and Table 0. We compare in § |^ the standard solar model with the helioseismological 
measurements and compare in § ^ the eight variant models with the observations . The 
dependence of the profiles of the sound speed, the density, and the adiabatic index upon 
the assumed reference model is determined in § Rl. We summarize and discuss our principal 
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2. Inversion technique 

The equations describing linear and adiabatic stellar oscillations are known to be 
Hermitian (cf. Chandrasekhar 1964). This property of the equations can be used to relate 
the differences between the structure of the Sun and the reference model to the differences 
in the frequencies of the Sun and the model by known kernels. Non-adiabatic effects and 
other errors in modeling the surface layers give rise to frequency shifts (Cox & Kidman 
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1984; Balmforth 1992) that are not accounted for by the variational principle. In the 
absence of a more fundamental method, these surface effects have been treated by the ad 
hoc procedure of including an arbitrary function of frequency in the variational formulation. 

When the oscillation equation is linearized — under the assumption of hydrostatic 
equilibrium — the fractional change in the frequency can be related to the fractional 
changes in the squared sound-speed (c^) and density (p). 

The sound-speed c used here is the adiabatic sound speed, which is defined as 



•2 



P 

where Fi, the adiabatic index, is defined by the thermodynamic relation 



(1) 



r. ^ ('^] . (2) 



c^lnp 



Here P is the pressure and s is the entropy. 
We can write, 



-= / A-k,(r)^(r)* + / K;,Ar)"Ar)dr + , (3) 

OJi J C J p Hj, 



(cf, Dziembowski et al. 1990). Here Sui is the difference in the frequency Ui of the ith 
mode between the solar data and a reference model. The kernels p and ^2 are known 
functions of the reference model which relate the changes in frequency to the changes in 
and p, respectively; and Ei is the inertia of the mode, normalized by the photospheric 
amplitude of the displacement. The term -Fgurf results from the near-surface errors. 

In this work we have used the subtractive optimally localized averages (SOLA) method 
(cf, Pijpers & Thompson 1992) to invert Eq. (|) in order to determine 5(? /(? and 5pl p 
between a reference model and the Sun. The principle of the inversion technique is to 
form linear combinations of Eq. (^ with weights di{ro) chosen so as to obtain an average 
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of 8(? !(? (or 5p/ p) localized near r = rg, while suppressing the contributions from 5p/p 
when inverting for 5c^/c^ (or 5(?/c^ when inverting for 6p/ p ) and the near-surface errors. 
In addition, the statistical errors in the combination must be constrained. We define the 
averaging kernel as 

/C(ro,r) = ^d,(ro)ir^.,,(r), (4) 

which is normalized so that / }C{rQ,r)dr = 1. The results of the inversion determine, for 
example, the difference, 6c^ between the square of the model sound speed and the square of 
the sound speed of the sun, i. e.. 

If the inversion is successful, the relative sound-speed difference can be written as 

/Sc^ Soj' 
K:{ro,r)—dr c::i^di{ro) — -. (6) 

The second expression is only approximately equal to the third expression in Eq. (^ because 
contributions from the second and third terms in Eq. (^ cannot be eliminated completely. 
The weights di are determined such that these contributions are substantially less than the 
error in the solution because of measurement errors in the data. An expression similar to 
Eq. (H) can be written for (^p/p);^^ ('"o) as well. 

The averaging kernels, /C, determine the extent to which we obtain a localized average 
of 5c^/(? or 5p/p. The width of the averaging kernel, e.g., the distance between the first 
and third quartile points, provides a measure of the resolution. Ideally, one would like 
the averaging kernel to be a 5-function at r = tq, but since only a finite amount of data 
is available, that is impossible to achieve. The effect of finite resolution on the inferred 
values of be? je? was studied by Bahcall, Basu & Kumar (1997) and found to be small for 
contemporary data sets. 
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The errors in the inversion results are calculated assuming that the errors in the 
frequencies are uncorrelated. Thus the error in (5c^/c^)inv at radius tq is given by 

a\r,)=j:dKro)-\ (7) 

where, is the quoted error of mode i with frequency uji. 

The adiabatic index Fi is related to the sound speed and density [cf., Eq. Hence, 
the kernels of and p can be easily converted to those of Fi and p (and vice versa, see e.g., 
Gough 1993). 

The details of how the method is implemented can be found in Basu et al. (1996) and 
the effect of the inversion parameters on the results are discussed in Rabello-Soares, Basu 
& Christensen-Dalsgaard (1999). 



3. Models used 

We have used a set of nine solar models as reference models to invert the three sets 
of helioseismological data. Figure |l| compares the computed sound speeds and densities of 
eight variant and deficient solar models discussed in § ISl^ and § O] with the sound speeds 



and densities computed for our standard model, STD, which is described in § 3.1 



Table |l| summarizes some of the key properties of the solar models discussed in this 
section. The convective zone (CZ) for most of the models are very close to the observed 
value of 0.713 ± 0.001 (cf, Basu & Antia, 1997). The only models for which the depth 
of the convective zone is clearly wrong are the NODIF and MIX models. The surface 
helium abundance of most of the models is also consistent with the abundance determined 
hehoseismologically, 0.248 ± 0.003 (Basu 1998). Only the NODIF model has a helium 
abundance that is obviously inconsistent with the observed helium abundance. 
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3.1. Standard Model: STD 

Our standard model, STD, is constructed with the OPAL equation of state (Rogers, 
Swenson & Iglesias 1996) and OPAL opacities (Iglesias & Rogers 1996) which are 
supplemented by the low temperature opacities of Alexander & Ferguson (1994). The 
model was constructed using the usual mixing length formalism for calculating convective 
flux. The nuclear reaction rates were calculated with the subroutine exportenergy.f (cf. 
Bahcall and Pinsonneault 1992), using the reaction data in Adelberger et al. (1998) and 
with electron and ion weak screening as indicated by recent calculations of Gruzinov & 
Bahcall (1998; see also Salpeter, 1954). The model incorporates helium and heavy element 
diffusion using the exportable diffusion subroutine of Thoule (cf. Thoule, Bahcall & Loeb, 
1994; Bahcall & Pinsonneault 1995)[[ 

For the standard model, the evolutionary calculations were started at the main-sequence 
stage. The model has a radius of 695.98 Mm. The ratio of heavy elements to hydrogen 
{Z/X) at the surface of the model is 0.0246, which is consistent with the value observed by 
Grevesse & Noels (1993). A Krishnaswamy T-r relationship for the atmosphere was used. 

Earlier inversions with similar models have shown that the difference in sound-speed 
between standard solar models and the Sun is small, of the order of 0.1% rms (cf. Bahcall 
et al. 1997; Basu et al. 1997; Bahcall, Basu & Pinsonneault 1998). 



Both the nuclear energy generation subroutine, exportenergy.f, and the diffusion 
subroutine, diffusion. f, are available at the Web site www.sns.ias.edu/ jnb/,menu item: 
neutrino software and data. 
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3.2. Variant models 

In this subsection, we describe four models that are slight variations on the theme of 
the standard model. 

Model PMS is evolved from the pre-main sequence stage, but otherwise is the same as 
STD. The difference in internal structure that results from including the pre-main sequence 
evolution is known to be very small (Bahcall & Glasner 1994). 

Model ROT incorporates mixing induced by rotation and is a reasonable upper bound 
to the degree of rotational mixing which is consistent with the observed depletion of 
lithium in the sun (cf. Pinsonneault 1996, Pinsonneault et al. 1999). The initial rotation 
period for the model is 8 days, which is the median observed value for T Tauri stars 
(Choi & Herbst 1996). The structural effects of rotation were treated using the method 
of Kippenhahn & Thomas (1970) Rigid rotation as a function of depth was enforced at 
all times in convective regions; in radiative regions the transport of angular momentum 
and the associated mixing were solved with coupled diffusion equations (see section 5 
in Pinsonneault 1997). Angular momentum loss from a magnetic wind is included, and 
then the thermal structure and angular momentum distribution in the interior are used to 
determine the diffusion coefficients as a function of depth. The angular momentum loss rate 
and the velocity estimates for the diffusion coefficients are the same as Krishnamurthi et al. 
(1997). The parameters of the model, which is started in the pre-main sequence phase, are 
fixed as follows: (1) the model is required to reproduce the equatorial surface rotational 
period of 25.4 days; and (2) the model is required to reproduce the observed solar lithium 
depletion of 2.19 dex (the difference between the meteoritic Li abundance of 3.34 on the 
logarithmic scale where H=12 and the photospheric Li abundance of 1.15). The rotational 
model neglects angular momentum transport by internal magnetic fields and gravity waves; 
both of these mechanisms can transport angular momentum without mixing and therefore 
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reduce the angular momentum content in the core, decreasing the mixing from meridional 
circulation and different instabilities. There is evidence from helioseismology that additional 
angular momentum transport mechanisms, such as gravity waves or magnetic fields, are 
needed to explain the absence of strong differential rotation with depth in the solar core 
(see Tomczyk, Schou & Thompson 1995). 

There has been considerable discussion recently regarding the precise value of the solar 
radius (cf. Antia 1998; Schou et al. 1997; Brown & Christensen-Dalsgaard 1998) and some 
discussion of the effects of the uncertainty in radius on the quantities inferred from the 
helioseismological inversions (cf. Basu 1998), We have therefore considered two models 
which were constructed with the same input physics as STD, but which have model radii 
which differ from the radius assumed in constructing STD. 

Model R78 has a radius of 695.78 Mm, which is the radius that has been determined 
from the frequencies of f-modes (cf. Antia 1998). Model R508 has a radius of 695.508 Mm, 
which is the solar radius as determined by Brown & Christensen-Dalsgaard (1998), who 
used the measured duration of solar meridian transits during the 6 years 1981-1987 and 
combined these measurements with models of the solar limb-darkening function to estimate 
the value of the solar radius. The solar structure is affected only very slightly by the choice 
of model radii. The fractional differences in the model radii considered in this paper are 
less than 1 part in 10^, whereas the radial resolution in the sound speed is at best a few 
percent (see Figure 1 of Bahcall, Basu & Kumar 19u97). 

The rms sound speed differences between the variant models and the STD model are: 
0.03% (PMS), 0.08% (ROT), 0.15% (R78), and 0.03% (R508). The average difference (rms) 
between the four variant models and the STD model is 0.07% . 
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3.3. Deficient models 

In this subsection, we describe four models that are each deficient in one or more 
significant aspects of the input physics. 

The model, OLD, is a standard solar model constructed with some relatively old 
physics: the Yale equation of state (cf. Guenther et al. 1992) with the Debye-Hiickel 
correction (cf. Bahcall, Bahcall & Shaviv 1968) and old OPAL opacities (Iglesias, Rogers 
& Wilson 1992) supplemented with low temperature opacities from Kurucz (1991). The 
model does include helium and heavy element diffusion and uses the nuclear reaction cross 
section factors (Sq) from Adelberger et al. (1998). In the course of writing this paper, we 
uncovered a small inconsistency in the code for the Yale equation of state. Fortunately, this 
inconsistency(which was introduced in recent revisions) does not affect any of our published 
results which no longer use the Yale equation of state. For the OLD model, the error in the 
code causes an increase in the mean molecular weight at fixed composition of 0.1% relative 
to the correct value. 

The OLD model differs from the STD model in using a cruder equation of state and less 
precise radiative opacities. Using the old physics rather than the current best input data as 
was done in constructing STD causes significant changes, primarily in the convection zone. 
This is a typical signature for large differences in the input equation of state, which is the 
most significant physics deficiency of this model. 

For Model S34 = 0, the cross-section of the nuclear reaction '^Hc(q;,7)^Bc was set equal 
to in order to minimize the calculated neutrino capture rates in the Gallex and SAGE 
experiments (see Bahcall 1989, Chapter 11). This assumption contradicts many laboratory 
experiments which measured a cross section for the ^He(Q;,7)'''Be that is competitive with 
the other way of terminating the p — p chain, namely, ^He(^He,2p)^He . For the S34 = 
model, nuclear fusion energy is achieved in a significantly different way than for the 
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standard solar model and therefore the calculated solar structure is appreciably different 
from the standard model (Bahcall and Ulrich 1988). In the standard solar model, about 
15% of the terminations of the p — p chain involve the '^He(a, 7) ''Be reaction, whose rate is 
proportional to S34. If we artificially choose S34 = 0, then in this model ''Be is not produced 
and there are no ''Be or ^B neutrinos. 

Model NODIF does not include either helium or heavy-element diffusion. This model 
therefore represents the state-of-the art in solar modeling prior to 1992 (cf. Bahcall & 
Uhich 1988, Bahcall & Pinsonneauh 1992, Proffitt 1994). 

Model MIX has an artificially mixed core. The inner 50% by mass (25% by radius) 
was required to be chemically homogeneous at all times. All of the other ingredients of the 
model, including helium and heavy element diffusion, are the same as in model STD. This 
model was constructed to be similar to the prescription of Gumming & Haxton (1996), who 
changed by hand the ^He abundance as a function of radius in the final BP95 (Bahcall & 
Pinsonneault 1995) solar model in order to minimize the discrepancy between measurements 
of the total event rates in neutrino experiments and the calculated event rates. Since the 
sun evolves over time and Gumming and Haxton only changed the abundances in the final 
model, we had to adopt some additional prescription as to how the mixing proceeds as a 
function of time. We assumed that the mixing was infinitely effective (the core was fully 
mixed) and constant in time. 

The rms sound speed differences between the deficient models and the STD model are: 
0.4% (NODIF), 0.2% (OLD), 0.2% (^34 = 0), and 1.9%(MIX). The average rms difference 
between the deficient models and the STD model is 0.7% , which is an order of magnitude 



larger than for the variant models discussed in § |3]2 . 
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4. How accurate are the measurements? 

How similar are the sound speeds, densities, and values of Fi inferred from different 
data sets? For sound speeds, this question is answered in Fig. |^(a)-(c), which shows the 
sound-speed difference between the standard model STD and the Sun as obtained using the 
MDI, GONG, and BiSON+LOWL data sets. The results appear relatively similar to the 
eye, but there are some differences as large as 0.05% inside 0.2Rq. For densities and Fi, the 
differences between the standard model and the Sun are illustrated in Fig. |^. 

4.1. The sound speed 

Figure ^ and Figure reveal even more clearly the differences between the three 
data sets. In Fig. | we plot the difference in sound speeds obtained with the GONG and 
the MDI data and the difference in sound speeds obtained with BiSON+LOWL and MDI. 
In all the panels shown in Fig. |^ and Fig. only one solar model, STD, has been used. 
The resolution of the inversions obtained with the three data sets is almost the same, hence 
the errors due to resolution should be very similar in each set. The only exception occurs 
near the solar surface, where the fact that the three data sets have different high-degree 
coverage becomes important. The extent of the high-degree coverage is probably the cause 
of the systematic differences in the sound speeds that are seen in Fig. ^ and Fig. ^ near the 
surface. Elsewhere, differences are caused solely by the measurements. 

The errors in the velocity measurements are apparently reasonably well understood; 
nearly all of the points lie within the 2a error bounds delineated in Fig. |[ If there were 
large systematic errors in one of the experiments, then we would have expected to see values 
for Ac/c in Fig. ^ that fell outside the 2a limits. The errors shown in Fig. |^ were calculated 
by combining quadratically at each target radius the errors obtained for inversions of the 
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MDI set and the other sets. The errors at each radius were evaluated as per Eq. 0. 

All data sets yield results for the sound speed profile that are consistent with each 
other within the errors of the data sets (cf. Fig. The rms differences are only 0.02% 
for the sound speeds calculated with the BiSON + LOWL and the MDI data sets and also 
0.02% for the differences found between the sound speeds calculated with the GONG and 
the MDI data sets. 

4.2. The density profile 

The density profile cannot be determined as precisely as the profile of the sound speed. 
The primary reason for the reduced precision in inverting for the density profile is that the 
p— mode oscillation frequencies are determined predominantly by the sound speed, with the 
density contributing only through the perturbation of the gravitational potential. In fact 
in the asymptotic limit, the frequencies are determined by sound-speed alone. As a result 
there is relatively little information about density in the frequencies. A further difficulty 
that must be overcome in a density inversion is the precise satisfaction of the condition for 
mass conservation, 

/ A7cr^6p = 0. (8) 
Jo 

In order to satisfy this condition with high numerical accuracy, the density must be 
reasonably well determined at all radii. Equation (|^) therefore requires a set of oscillation 
frequencies that includes a large number of low degree modes (to invert accurately for the 
core ) as well as a large number of high degree modes (to invert accurately for the surface). 
If a proper set of either high degree or low degree modes is not available, the density 
inversion becomes very uncertain. 

Only for the MDI data were we successful in forming a local averaging kernel that 



-17- 



permitted a good inversion for solar density. Fig. ^(a) shows that the difference between the 
density profiles of the solar models and the helioseismologically determined density profile 
is < 1% . However, the accuracy of this measurement is an order of magnitude less precise 
than for sound speeds. 

4.3. The adiabatic index Fi 

The Fi difference obtained between model STD and the Sun is shown in Fig. ^(b) . 
Since most of the Fi difference between the models and the Sun is concentrated at the 
surface, we use only the MDI data for the Fi inversions. This set has the most number of 
high degree modes. 

With the convenient inversion method used here (SOLA), we do not have good spatial 
resolution close to the surface. There are computationally intensive inversion methods 
which give superior resolution near the surface, for example. Optimally Localised Averages 
(cf. Kosovichev 1995, Elliot & Kosovichev 1998; Basu, Dappen & Nayfonov 1999). For the 
general survey performed in this paper, we did not feel it necessary to invert very close to 
the solar surface. 

The contribution to the solution for Fi from the second function (density in this case) 
is more difficult to suppress than for inversions of sound speed or density. Thus we expect 
larger errors for reference models with large density differences relative to the Sun. 

Given the greater precision of the sound-speed measurements, we will emphasize in 
what follows the profile of the sound speed and will only refer to the density and Fi profiles 
for completeness. 
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5. How large are the effects of finite radial resolution? 

We calculate the solar sound speed using the relative sound-speed difference between 
the models and the Sun. Thus, if {5c^/c'^)mv is the result of the inversion, then 



C0)inv(ro) 



(r„) + l 

V / mv 



Cmodel(^o)- (9) 



However, {cQ)iay{ro) is not identically equal to the true solar sound-speed, since the 
inverted speed is an average of the sound-speed difference in the vicinity of Tq. The 
averaging kernel at Tq, /C(ro, r) defines the region in r over which the averaging is done. We 
want to estimate the error introduced in the inferred sound speed due to finite resolution. 
Since we do not know the true sound speed inside the Sun, we estimate the errors using a 
solar model as a "proxy Sun." Thus if Cproxy is the sound speed of the proxy Sun, Cmodei is 
the speed in the reference model, and 5(? /(? = (Cp^oxy — Cmodci)/ 'Cproxy the relative difference 
between the squares of sound speeds of the two solar models , then the relative error in the 
inferred sound-speed of the proxy Sun due to the finite resolution of the averaging kernel is 



5c\ . \ 5c' 



(10) 



The factor of (1/2) in the above expression arises from the conversion of relative errors in 
(? to relative errors in c. Note that (Ac/c(ro))resoi = if the averaging kernel /C(ro, r) is a 
delta function. The error in density due to the resolution of the density kernel is similarly 
given by 

f^(ro)) =f/Wr„,r)^W*)-^(r„). (11) 

V P /rcsol V'' P ) P 

where /Cden is the averaging kernel obtained for density inversions. 

In the subsequent discussions, the profile of the solar sound-speed obtained with MDI 
data using STD as the reference model is referred to as the "standard sound speed profile" . 
Similarly the solar density profile inferred from MDI data using STD is referred to as the 
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"standard density profile," and the solar Fi profile inferred MDI data using STD is referred 
to as the "standard Fi profile." 

Figure |^ shows (Ac/c)resoi for model NODIF and ROT when STD is used as the proxy 
sun. The resolution errors are small and generally less than 0.02-0.03 percent. However, 
the resolution errors are relatively large in the solar core and the base of the convection 
zone. The large error in the core results from the fact that there are very few p-modes 
that sample this region well. This causes the averaging kernel to be relatively wide in this 
innermost volume. The even larger resolution error at the base of the convection zone is 
caused by the sharp gradient in the 5(? /(? for the different models. 

The rms difference between the sound speed profile inferred with the ROT model as the 
reference model and the STD sound profile as the proxy sun is 0.015%; the rms difference 
is 0.038% when the NODIF model is used as the reference model. The rms difference for 
the density profile is 0.080% when ROT is the reference model and 0.31% when NODIF is 
used as the reference model. 

The errors are larger when NODIF is compared with STD than when ROT is compared 
with the STD model, which simply reflects the fact that the difference in sound speed 
profiles between model STD and NODIF is larger than the difference in sound speed proflles 
between STD and ROT. This observation leads to the rather obvious conclusion that one 
expects to get more accurate solar sound speeds by using reference models which have 
sound speeds that are similar to the Sun. 

The errors in the inferred density are also shown in Fig. ^. The errors in density are 
about an order of magnitude larger than the errors in the sound speed. 

The errors in the inferred Fi due to resolution effects are expected to be small. The 
reason is that in the region we have been able to do the inversions (r < 0.94 R©), Fi 
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differences between most models and the Sun are very smooth (see for example Basu et al. 
1999) The only model which could cause large errors is model OLD. 

6. Standard sound speeds, densities and Fi 

Table ^ lists the solar sound speed, density, and Fi profiles that are obtained with the 
STD model using MDI data. These data may be useful for other applications. Therefore, 
we have made available machine readable files in the format of Table |^, but with a denser 
grid in radius, at the web-site [http:/ /www. sns.ias.edu/~jnb . 

7. Comparison with the standard solar model 

How well does the standard solar model, STD, agree with the different measurements 
of the sound speed? Is the difference, shown in Fig. ^(a)-(c), between the STD model and 
each of the measurements larger or smaller than the differences between the measurements 
themselves (shown in Fig. ^)? 

The root-mean-squared difference between the STD sound speeds and the solar speeds 
is 0.069% for the MDI data, 0.069% for the GONG data, and 0.064% for the BiSON+LOWL 
data set. These results are averaged over all regions of the sun for which good data are 
available, from O.O5-R0 to 0.95Rq. In all cases, the agreement is excellent, although the STD 
model can and should be improved, especially near the base of the convective zone (see, 
e.g., the discussion below of the pre-main sequence model, PMS, and the model including 
rotation, ROT). In the solar core, where the neutrinos are produced, the rms agreement is 
even slightly better: 0.062% for the MDI data, 0.061% for the GONG data, and 0.044% for 
the BiSON+LOWL data set. 
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For all three data sets, the sound speeds at the base of the convection zone of the STD 
model differ by about 0.2% from the helioseismological values . This discrepancy has been 
seen earlier with similar standard models (cf. Gough et al. 1996, Bahcall et al. 1997) and 
can be attributed to the lack of mixing in the model below the base of the convective zone. 
In the models, diffusion without mixing causes a sudden, local rise in the helium abundance 
below the base of the convective zone. The increase in helium abundance increases the 
mean molecular weight, thereby decreasing the sound speed (since (? oc T//i, where T is the 
local value of the temperature and [i is the local mean molecular weight). 

The rms density difference between the model and the Sun is within 1% [see Fig. ^a)]. 
For Fi, the rms difference between the STD model and the Sun is less than 0.1% [Fig. 0(b)] . 
The somewhat larger difference in the core has been attributed to errors in the equation of 
state (cf., Elliott & Kosovichev 1998). 



8. Comparison with variant and deficient models 



How well do the variant models discussed in § |3.2| and § p.3| agree with the 
helioseismological measurements? 

Figures ^ and show the results of the inversions made using, respectively, MDI 
(Fig. §) and GONG and BiSON+LOWL (Fig. |7|) data. For most of the models, the vertical 
scale for the fractional velocity differences has a range of a few tenths of a percent. For the 
densities, the corresponding range is of order a few percent. 

For Fi, Fig. ^ shows the fractional differences between the Sun and the solar models 
that were found using MDI data. The fractional differences are less than or of order of 0.2% 
for all of the solar models except for the MIX model. For the MIX model, the fractional 
differences are larger, of order 0.5% . We conclude that, with the exception of the MIX 



-22- 



model, the theoretical profiles of Fi are in good agreement with the solar values of Fi. 

Table ^ summarizes the average root-mean-squared deviation between the predicted 
sound speed profile of the solar models discussed in § ^ and the measured sound profile 
determined with the MDI data. 

8.1. Variant models 

The pre-main sequence PMS model yields values for 5c/ c that are similar to the results 
obtained with the STD model (comparing Fig. ^ and Fig. |^a,e with Fig. |^-c). This 
similarity is to be expected since the difference between the STD and PMS models is small, 
of order hundredths of a percent in 5c/ c everywhere and of order a few tenths of a percent 
in 5p/ p (cf. Fig. 0a,d). The difference in Fi between the PMS model and the Sun is similar 
to that between the STD and the Sun [cf. Fig. ^b) and Fig. Ka)]. 

The model with rotational mixing, ROT, agrees better than the STD and PMS 
models with the helioseismological measurements near the base of the convective zone (see 
Figures Ha,e and 0a,e). This improved agreement confirms the suggestion (see Richard et 
al. 1996) that mixing at the base of the convection zone is a possible explanation for the 
significant discrepancy in this region between the measured and the STD model sound 
speeds. However, with this version of mixing, the agreement is slightly worse in the solar 
core, resulting in an overall rms deviation that is essentially the same as for the STD model. 
The difference in Fi between ROT and the Sun is very similar to differences in Fi found 
with models STD and PMS. 

The two models with slightly different radii, R78 and R508, yield results (see Table ^ 
for the rms agreement with the MDI data that are comparable to the STD model. The R78 
model yields slightly better agreement and the R508 model yields slightly worse agreement 
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than is obtained with the standard model radius. The shape of the sound speed differences 
between the R78 and R508 models and the Sun is very similar to the shape of sound speed 
differences between the STD model and the Sun, but for the 'non-standard' radii the sound 
speeds are shifted downward in Fig. ^jb and Fig. ^,f (see also Fig. |c,d). 

8.2. Comparison with deficient models 

The model constructed using old input data, OLD, produces a significantly worse rms 
sound speed discrepancy, 0.17%, compared to the 0.07% for the STD model(see Table ^. 
It is encouraging that the improvements in nuclear physics, equation of state, and radiative 
opacity that are described in § have resulted in better agreement, by about a factor 
of two, with the measured sound speeds. The OLD model also shows larger difference, 
relative to the Sun, in Fi towards the surface, confirming our suspicion that the equation 
of state used in the OLD model is not sufficiently accurate for optimal helioseismological 
applications. 

The model with S34 = does not allow the nuclear reaction '^He(a,7) ''Be . This 
change in the nuclear physics results in a sufficiently large modification in the structure of 
the solar core in the model that the difference is easily seen in precise helioseismological 
measurements (see Bahcall et al. 1997). Figures |c and |^,g show that the sound speed 
predicted by the S34 = model differs from the helioseismologically inferred sound speeds 
by as much as 0.5% in the solar core, about an order of magnitude worse agreement than is 
obtained for the core with the STD model. The influence of the reaction is seen even more 
dramatically in the density differences with respect to the Sun. Figure ^(g) shows that the 
discrepancies in density are as large as 5% in the outer region of the sun. This result can 
be understood as follows. To achieve the same luminosity, the density in the core must be 
increased when an important nuclear reaction, ^He(a,7) ''Be, is artificially set equal to zero. 
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Since mass is conserved, any change in the density in the core has to be compensated by 
an opposite, larger change in the less dense outer layers. The small reduction relative to 
the STD model in the density in the core of the S34 = solar model results in a relatively 
large change in the density in the outer layers and a significant discrepancy with the 
helioseismologically inferred density profile. Since the equation of state used in this model is 
the same as in STD, the differences in Fi between this model and the Sun are very similar 
to those found with the STD model. 

The model NODIF is a generally poor fit to the helioseismological measurements. 
Figures ^,h and |^,h show that the disagreement is consistently large near the base of 
the convective zone, which refiects the fact that omitting diffusion results in models with 
incorrect depth of the convection zone (Bahcall & Pinsonneault 1995, cf. also Table |^ of the 
present paper). The rms discrepancy between the sound speeds of the NODIF model and 
the measured sound speeds is 6.5 times worse on averaged over the sun than for the STD 
model. In the region where we have inverted for Fi difference, NODIF fares quite well. The 
main difference in Fi between the NODIF model and the Sun is expected to arise because 
of differences in helium abundance. But that difference will show up only in the helium 
ionization zone (around O.98i?0) which we have not resolved. 

The vertical scale for 5c/c must be increased by a factor of 25, from 0.002 to 
0.05(negative discrepancy), in order to display the very large discrepancy that exists for the 
MIX models. The differences are particularly glaring in the core, where the model is fully 
mixed. The sound speed difference between the MIX model and the measurements is as 
much as 5% in the solar core while the density difference is almost 40%. The corresponding 
maximal differences for the standard model in the solar core are 0.1% in 5c/c and 1% in 
(5p/p, about 50 times and 40 times smaller than for the MIX model. The average rms 
discrepancy for the MIX model is about 25 times larger than for the STD model. Obviously, 
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this model is not a good model of the sun, although it was proposed (Gumming & Haxton 
1996) as a way of decreasing (but not eliminating) the differences between standard neutrino 
flux predictions and the measured neutrino fluxes. We were unable to obtain a reliable 
inversion for Fi with model MIX because of the large difference in density between the 
model and the Sun. 



9. Dependence upon the assumed reference model 

In this section, we evaluate the dependence of the inferred sound speed profile and the 
inferred density profile upon the assumed reference model. We calculate the sound speed 
or the density using two different combinations of reference model and measurement data 
and then compare the results. Thus we evaluate the set of differences formed by (model, - 
datafc) - (modelj - dataj), where for convenience we always take model^ as the STD model 
and data-,- as the MDI data set, but model, can be any one of the eight variant models 
discussed in § O and data^ is either the MDI, BiSON+LOWL, or GONG data set. 



The principal results of this section are summarized in Table 



9.1. Dependence of sound speed profile upon reference model 



Figures ^ and |T0| show the relative differences (double differences in the sense defined 
above) between the standard sound speed profile (obtained with the model STD and MDI 
data ) and the solar sound speed inferred using the eight variant models. We present result 
for all three of the data sets: MDI, BiSON+LOWL, and GONG. 



The fractional difference, Ac/c, that is shown in Fig. ^ and Fig. |TO]is defined explicitly 



- 26 - 



by the relation 

Ac _ Cq^ variant ~ Cq, STD 
C Cq, STD 

The quantities Cq^ variant and Cq^ std are the best estimates for the solar sound speed 
that are found using the specified reference model and data set. More explicitly, 
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C0, model = Cmodoi\/(1 + / c'^)modeh whcrc {S / c^) model IS obtained dircctly from the 
inversion equation, Eq. (||). 

The inversions obtained using the PMS, the ROT, the OLD, and the S34 = models 
all yield sound speeds that differ from the standard sound speed profile by only a few 
hundredths of a percent, i. e., within the errors due to the measurements (see Table ^ and 
Fig. ^ and Fig. [T^). The only conspicuous exception to this statement occurs near the base 



of the convective zone, where the finite resolution causes a difference that is larger than the 
recognized measurement errors (see Fig. ^b). 

The largest systematic, monotonic differences, ~ 0.03% to 0.08%, are found for the 
models R78 and R508(see Fig. §c,d), although other models have larger rms differences (cf. 
Table ||). The sound-speed profiles obtained using these two models as reference models 
(and the MDI data for the measurements) exhibit a smooth difference with respect to the 
standard sound-speed profile. Data errors and resolution errors are not important when the 
MDI data are used for both the standard and the variant model inversions. One can see the 
irregular effects of using different data sets with finite resolution in the BiSON+LOW and 
GONG panels of Fig. Pc,d. Equation (0) is obtained assuming that there is no difference in 
radius between the Sun and the model; hence, a difference in the radius between the model 
and the sun can introduce a systematic error in the inversion results. One must use an 
accurate value for the solar radius in order to obtain precisely correct inversion results (see 
Antia 1998 and Basu 1998). 

Even when the NODIF model is used as the reference model, the inferred sound speed 
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profile is in reasonable agreement with the standard profile except near the base of the 
convective zone. At the base of the convective zone, the combined effect of finite resolution 
and non-linear effects increase the velocity difference well beyond what is expected from 
measurement errors alone. The results for the NODIF model are shown in Fig. |lO|c. 

The linear inversion fails for the MIX model in the solar core and near the base of 
the convective zone (see Fig. p!0|d). The rms fractional differences, Ac/c, are ~ 0.4% (see 
Table an order of magnitude larger for the MIX model than for the models which 
more closely resemble the STD model. This failure is not surprising since the MIX model 
is very different from the sun (and the STD model) in the core and at the base of the 
convective zone. What is more remarkable is that despite the very large difference between 
the reference model and the Sun, the linear inversion scheme yields sound speeds for much 
of the solar volume that are within a few tenths of a percent of the results obtained using 
much better reference models. 



9.2. Dependence of the inferred density profile upon reference model 

Figure |TT| shows the dependence of the inferred density profile upon the assumed 
reference model. We are able to make these comparisons only for one data set, the MDI 
data set, since we were not able to make satisfactory inversions for the density profile using 
the other data sets. 

For the PMS, ROT, R78, R508, and OLD models, which differ from each other only by 
modest amounts, the dependence of the inferred density profile upon the assumed reference 
model is moderately large, of order 0.2% (see the next-to-the-last column of Table H), but 
is nevertheless generally smaller than the estimated measurement uncertainties. Of course, 
this dependence upon reference model is about an order of magnitude larger than the 
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dependence of the sound speeds upon reference models (see above). 



For the NODIF, S34 = 0, and MIX models, the dependences shown in Fig. |rT]are much 
larger than the measurement errors, i.e. they are ~ a few percent. Non-linear effects are 
clearly important in these inversions. Nevertheless, the most remarkable fact may be that, 
despite the very significant differences between the variant models (NODIF, S34 = 0, and 
MIX) and the STD model, the different models all yield estimates for the solar density that 
agree with each other within a few percent. 



9.3. Dependence of the inferred Fi profile upon reference model 



Figure 12 shows the dependence of the inferred Fi profile upon the assumed reference 



model. All the comparisons are for the MDI data set. 

For models PMS, ROT, R78, R508 and S34 = 0, the Fi profiles obtained agree well 
with the Sun and with each other ( ~ of a few hundredths of a percent, see the last column 
of Table 4 for the rms dependence upon the reference model). Although the structure of the 
S34 = model is quite different from the Sun (and model STD), the structural difference is 
not large enough to cause major problems with the Fi inversion. The Fi profile obtained 
with model NODIF shows a moderately large difference, 0.05%. The profile obtained with 
model OLD shows an even larger difference, which we believe is due to an inadequate 
description of the structure close to the solar surface. By far the largest differences are 
found for model MIX, which is due to the large difference in structure between this model 
and the Sun. 
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10. Discussion 

The principal purpose of this paper is to explore some of the systematic uncertainties 
that affect the determination of the profiles of the solar sound speed and the solar density. 

As a by-product of this investigation, we have confirmed that standard solar models 
are in remarkable agreement with helioseismological measurements of the sun. For example, 
the rms difference between the standard solar model profile for sound speeds and the 
helioseismological profile is only 0.07% (see discussion in § Including pre- main sequence 
evolution or a small amount of rotationally induced mixing does not affect the average 
results very much, but can give better agreement with observations near the base of the 
convection zone(see discussion in § |^. 

Table ^ shows that five state-of-the-art solar models(STD,PMS,ROT,R78,R508), each 
constructed with some different physics or parameter choice, all give comparable agreement 
with the global helioseismological measurements. On the basis of the global seismological 
evidence, one cannot say that one of these models is definitely more like the sun than the 
other models. 

10.1. Systematic uncertainties 

10.1.1. Uncertainties in the data 

We determine in § ^ the systematic differences due to the choice of the individual data 
set by comparing sound speed profiles calculated using different data sets. The results 
are shown in Fig. Q the difference between the results from state-of-the-art data sets is 
rms about 0.02% averaged over the sun and bounces around within the 2a error envelopes 
determined by the combined measurement errors. 
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10.1.2. Effects of finite resolution 

We estimate in § ^ the uncertainties due to the finite resolution of the inversion kernel 
by adopting a particular solar model as a proxy sun and then comparing the convolved and 
inverted sound speed (or density) profile with the true profile in the proxy sun. The finite 
resolution of the inversion kernel leads to rms systematic uncertainties in the range ~ 0.02% 
to 0.04% in the profile of the sound speed, although the errors are typically much larger in 
the solar core and at the base of the convective zone( see Fig. The uncertainties in the 
density profile due to finite resolution are typically an order of magnitude larger than the 
errors in the profile of the sound speed. 



10.1.3. Uncertainties due to reference models 

We use nine different solar models in order to determine the effects of the choice of 
reference model upon the inferred sound speed, density and Fi profiles (see the discussion 
in sect H). The results are summarized in Table ^ 

We have performed calculations for a standard solar model (STD) and four variant 
models(PMS,ROT,R78, and R508, each described in § |3.3|) . All five of these models include 
physics and input parameters that are at the state-of-the-art for 1999 solar models. The 
average rms difference between the sound velocities of each of the variant models and the 



STD model is 0.07% (see § |3.2| ). The average rms difference between the sound speed profile 
inferred for the Sun using one of the variant models and the STD model is 0.03%(averaging 
the first four lines of Table 1). Hence, the spread among the inferred solar sound speeds is 
more than a factor of two less than the spread among the reference models themselves. 

We also performed calculations for four deficient models(OLD, S'34 = 0, NODIF, and 
MIX). The physics used in constructing each of these models is deficient in some significant 
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way (see § 3^). These deficiencies are reflected in the fact that the average rms difference 
between the sound velocities of each of the four non-standard models and the STD model 
is 0.7% (cf. § p.3|) , an order of magnitude larger than for the variant models. Nevertheless, 
these deflcient-by- design models give, when used as reference models, reasonably accurate 
values for the inferred solar sound speeds. The average rms difference between the sound 
speed proflle inferred for the Sun using one of the deflcient models and the STD model is 
0.13% (averaging the last four rows of Table Thus the discrepancy, when averaged over 
the different deflcient models, is a factor of more than flve less than the spread among the 
reference models. 

Our bottom line on the systematic uncertainties for sound speeds is that, as expected, 
even relatively crude reference models yield reasonably good estimates for the solar sound 
speed. 

Table ^ and Table ^ show that the proflle of Fi is determined with a precision that is 
similar to, or slightly better than, the proflle of the solar sound speed, that is to an accuracy 
~ 0.1%. The density is determined with an order of magnitude less precision, ~ 1%. 



10.2. Can helioseismology rule out some non-standard solar models? 

What models are strongly disfavored (ruled out) by the helioseismological data? 
We choose as a flgure of merit (crudely analagous perhaps to 1 standard deviation) the 
largest rms difference in sound speed proflle between the Sun and one of the variant 
(state-of-the-art) solar models. This rms difference is 0.1% (for the R508 model, see 
Table The S34 = model has approximately twice as large a deviation (0.2%) and 
is therefore somewhat disfavored, but the OLD model (1995 physics) is perhaps still 
within the range of acceptability. Two models are strongly disfavored (ruled out at a high 
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significance level). The no diffusion model, NODIF, has a 0.45% rms difference in sound 
speed profilewith respect to the Sun; this is 4.5 times worse than the least successful of the 
variant models. The model with a mixed solar core, MIX(Cumming and Haxton 1996), has 
a rms difference of 1.8%, 18 times worse than the least successful of the variant models. We 
therefore conclude that the no diffusion model and the mixed model are ruled out at a high 
level of significance. 
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Fig. 1. — The relative sound-speed differences, 5c/ c, and the relative density differences, 
Sp/p, between the standard model STD and the other solar models. For each model, 
Sc/c= (csTD -Cmodci)/cmodci and 6p/p = (psTD -PmodoO/Pmodei- The models are described in 
Table |l]. The vertical scales for panels c), d), g), and h) are larger than for the other panels. 
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Fig. 2. — The relative sound-speed difference between the Sun and model STD obtained 
using three different solar oscillation data-sets [Panels (a),(b) and (c)]. The differences are in 
the sense (Sun — Model)/Model. The vertical error-bars indicate la errors in the inversion 
results because of errors in the data. The horizontal error-bars are a measure of the resolution 
of the inversion. 




Fig. 3. — The relative density differences [Panel (a)] and Fi differences [Panel (b)] between 
the Sun and model STD obtained using MDI data. The differences are in the sense (Sun — 
Model)/Model. 
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Fig. 4. — The relative errors in the solar sound speed obtained with model STD as the 
reference model using different data sets. The reference sound-speed is taken to be the one 
obtained by inverting MDI data. The difference with those obtained using BiSON+LOWL 
and GONG data were calculated. The 1 and 2a error envelopes due to errors in the data are 
shown as dotted lines. Since the reference model was the same for all three inversions, and 
since the resolution of the inversions using the three data sets are very similar, we do not 
expect any additional error due to finite resolution or to differences in the reference model. 
The rms differences are only 0.016% for the sound speeds calculated with the BiSON + 
LOWL and the MDI data sets and also 0.020% for the differences found between the sound 
speeds calculated with the GONG and the MDI data sets. 
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Fig. 5. — The error in the estimated sound-speed and density of two models due to tlie 
finite resolution of the inversions. In both cases the reference models used is STD. The 
quantities (Ac/c)resoi and (Ap/p)i.esoi, defined in Eq. |^ and[Tl| respectively, are estimates of 
the errors due to resolution that are expected in the solar sound-speeds and density obtained 
by inversion. The errors are largest where the sound-speed difference between the reference 
model and the test model shows a sharp gradient, as occurs for example, at the base of the 
convection zone. 
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Fig. 6. — The relative sound-speed and density differences between tlie Sun and different solar 
models obtained using MDI data. The differences are in the sense (Sun — Model) /Model. 
Horizontal error-bars arc not shown for the sake of clarity. The differences are larger for the 
models described in the bottom two sets of panels and therefore the vertical scales cover a 
wider range for these panels. The different models are described in Table 1. 
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Fig. 7. — The relative sound-speed differences between the Sun and different solar models 
obtained using GONG and BiSON+LOWL data. The differences are larger for the models 
described in the bottom two sets of panels and therefore the vertical scales cover a wider 
range for these panels. The different models are described in Table 1. 
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Fig. 8. — The relative Fi differences between the Sun and different solar models obtained 
using MDI data. This figure is similar to Fig, |] except that Fig. ^ refers to Fi instead of 
sound-speed and density. The vertical scale of the last two panels of Fig. ^ is larger than 
the other panels in order to accommodate the large difference for model MIX [Panel (h)] 
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Fig. 9. — The fractional difference between the solar sound speed inferred using the STD 
solar model as a reference model and the sound speed obtained using each of four variant 
solar models as a reference model. For the reference sound profile, the MDI data set was 
inverted. Differences in the inferred sound speeds are presented for all three data sets, MDI, 
BiSON+LOWL, and GONG, and for all of the variant solar models. The la error envelope 
due to data measurement errors is shown as the dotted line. 
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Fig. 10. — The same as Fig. but for four solar models with considerably different 
physics than the STD model. The fractional differences obtained using the MIX model 
as the standard model are an order of magnitude larger than for the other cases considered; 
therefore, the vertical scale for the bottom panel covers an order of magnitude larger range. 
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Fig. 11. — The relative errors in the solar density obtained using different reference models 
and the MDI data set. The reference solar density used to compute the differences is the 
one obtained using model STD as the reference model. The la error envelope due to data 
errors is shown as the dotted hne. The differences are larger for the models described in 
the bottom two sets of panels and therefore the vertical scales cover a wider range for these 
panels. 
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Fig. 12. — The relative errors in the adiabatic index Fi for the Sun. obtained using different 
reference models and the MDl data set. The reference solar density used to compute the 
differences is the one obtained using model STD as the reference model. The Icr error 
envelope due to data errors is shown as the dotted line. The differences are larger for the 
model described in the last panel (h) and therefore the vertical scales cover a wider range 
for panels (g) and (h). 
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Table 1: Properties of the solar models used. The first five models fisted all use input 
physics that is within the currently acceptable range. The last four models are all deficient 
in one or more aspects of the important input physics. 



Model 


T 

c 


Pc 


Ys 


^cz/ Rq 


Comments 




10^ K 


g cm~^ 








STD 


15.74 


152.98 


0.2453 


0.7123 


Standard model, incorporates diffusion, 












only main sequence evolution 


PMS 


15.72 


152.73 


0.2455 


0.7127 


Same as STD, but with pre-main sequence evolution 


ROT 


15.69 


148.30 


0.2530 


0.7155 


Same as STD, but with rotational mixing of elements 


R78 


15.73 


152.97 


0.2454 


0.7122 


Same as STD but with radius of 695.78 Mm 


R508 


15.73 


152.97 


0.2454 


0.7121 


Same as STD but with radius of 695.508 Mm 


NODIF 


15.44 


148.35 


0.2653 


0.7261 


Same as STD but with no diffusion 


OLD 


15.80 


154.52 


0.2470 


0.7111 


Old physics (see text) 


834 = 


15.62 


153.50 


0.2422 


0.7151 


Same as STD but with reaction constant S34 set to 


MIX 


15.19 


90.68 


0.2535 


0.7314 


Core mixing (see text) 
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Table 2: Solar sound speed, density, and adiabatic index, Fi, as derived from MDI data with 
model STD 
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Table 3: The root-mean-squared (RMS) 
sound speed, density, and adiabatic index 
differences with respect to the Sun. 
The table lists the percentage differences 
between the variables determined with the 
MDI data and those predicted by different 
solar models considered in the text (see 



Section 


3.2). The 


first four rows refer 


to variant standard 


models; the last four 


rows refer to different deficient 


(or non- 


standard) models. 


This table shows how 


different 


each of the models is 


from the 


Sun. 








Model 


(5c/c) 


(5P/P) 






% 


% 


% 


STD 


0.069 


0.942 


0.087 


PMS 


0.085 


0.929 


0.082 


ROT 


0.069 


1.404 


0.084 


R78 


0.064 


0.673 


0.090 


R508 


0.098 


0.408 


0.126 


OLD 


0.170 


0.905 


0.126 


834 = 


0.209 


3.108 


0.146 


NODIF 


0.447 


5.303 


0.093 


MIX 


1.795 


13.319 


0.341 
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Table 4: Dependence upon reference model: This table 
gives the rms differences in solar sound-speed, density, 
and adiabatic index obtained using different reference 
models between radii of 0.07Rq and 0.9Rq. See the 
introduction to Sect. ^ for the definition of how the 
reference model dependence is calculated. 



Model 


(^) 


















^ P ' 






MDI 


BiSON+LOWL 


GONG 


MDI 


MDI 




% 


% 


% 


% 


% 


PMS 


0.0246 


0.0256 


0.0131 


0.2245 


0.0124 


ROT 


0.0179 


0.0321 


0.0373 


0.3067 


0.0251 


R78 


0.0253 


0.0345 


0.0415 


0.1361 


0.0150 


R508 


0.0685 


0.0700 


0.0781 


0.3077 


0.0476 


OLD 


0.0190 


0.0250 


0.0125 


0.2206 


0.1401 


S34 = 


0.0224 


0.0251 


0.0226 


1.5556 


0.0231 


NODIF 


0.0734 


0.0465 


0.0826 


1.8858 


0.0543 


MIX 


0.4112 


0.5078 


0.3948 


3.3593 


0.2899 



